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Abstract 

In a recent article the most general non-uniform reaction-diffusion models 
on a one-dimensional lattice with boundaries were considered, for which 
the time evolution equations of correlation functions are closed and the 
stationary profile can be obtained using a transfer-matrix method. Here 
models are investigated for which also the equation of relaxation towards 
the stationary profile could be solved through a similar transfer-matrix 
method. A classification is given, and dynamical phase transitions are 
studied. 
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1 Introduction 



Most studies on reaction-diffusion models are focused on uniform lattices where 
reaction rates are site-independent. Among the simplest generalizations beyond 
a uniform system is a lattice with alternating rates. In [TH3], the steady state 
configurational probabilities of an Ising spin chain driven out of equilibrium 
by a coupling to two heat baths has been investigated. An example is a one- 
dimensional Ising model on a ring, in which the evolution is according to a 
generalization of Glauber rates, such that spins at even (odd) lattice sites expe- 
rience a temperature {To). In [4j the dynamics of a reaction-diffusion model 
with spatially alternating rates has been discussed. The response function to an 
infinitesimal magnetic field for the IsingGlauber model with arbitrary exchange 
couplings has been studied in In [5^, relaxation in the kinetic Ising model 
on an alternating isotopic chain has been discussed. The critical dynamics of 
a disordered Ising ferromagnetic chain with two coupling constants for Glauber 
dynamics was studied in [7]. 

In [8], a non-uniform extension of the Glauber model on a one-dimensional 
lattice with boundaries was investigated. In this model the reaction rates are 
based on detailed balance. There the expectation values of spins in an Ising 
model with nonuniform coupling constants were studied. A transfer matrix 
method was used to study the steady state behavior of the system in the ther- 
modynamic limit. Different (static) phases of this system were studied, and a 
closed form was obtained for this transfer matrix. In [9 , the dynamics and pos- 
sible dynamical phase transitions of the non-uniform extension of the Glauber 
model were investigated. It was shown that there are two phases: the fast phase 
and the slow phase. In the fast phase the relaxation time is independent of the 
reaction rates at the boundaries. Changing continuously the reaction rates, the 
system may experience a phase transition. In the slow phase the relaxation time 
changes with reaction rates at boundaries. 

In [lOj the most general non-uniform exclusion nearest neighbor reaction- 
diffusion models on a one-dimensional lattice with boundaries, were considered, 
for which the evolution equations of the one-point functions are closed and the 
stationary profile of the system can be obtained using a transfer-matrix method. 
It was shown that the criterion that the evolution equations be closed is the same 
as the case of uniform models. Such models are called autonomous. To be able 
to obtain a closed form for the transfer matrix, however, additional conditions 
are to be satisfied. Models satisfying these additional conditions as well, are 
called superautonomous. 

In this article we want to study the dynamics and possible dynamical phase 
transitions of superautonomous models, using a transfer-matrix method. The 
scheme of the paper is as follows. In section 2, a brief review of the formal- 
ism is given. In section 3, the evolution equations governing the expectation 
values of the number operators are studied. It is seen that solvability through 
the transfer-matrix method requires additional criteria, apart from those of 
the superautonomy conditions obtained in [TU]. In section 5 a classification is 
introduced for superautonomous models for which the relaxation towards the 
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stationary profile of tlic of expectation values of the number operators is solvable 
through the transfer-matrix method. In section 5, some examples are explicitly 
discussed. Section 6 is devoted to the concluding remarks. 



2 Exclusion nearest- neighbor reaction-diffusion 
models with nonuniform reaction rates 

Consider a one-dimensional lattice with {L + 1) sites, numbered from to L. 

Each site is either empty or occupied with one particle. The evolution of the 
system is governed by a Hamiltonian "H of the form, 

n = K+(^nc^+n'r^, (i) 

where "Ha corresponds to the link a: 

•Hg = 1 (g) ■ ■ -^(g) ■ ■ ■ 1 Ig) gg Ig) 1 (g) ■ ■ -^(8) ■ ■ ■ 1 , (2) 

a—fj. L — a—fj, 

and 

,:=l. (3) 

The link a links the sites (a — fi) and (a + /x), so that a ± /x are integers. Ha 
is a linear operator acting on a four dimensional space (the configuration space 
corresponding to the sites (a — /i) and (a + /x)) with a basis {eoo, eoi, eio, ei i}, 
where denotes vacancy and 1 denotes occupied site. Also, 

n'o = Hq (g 1 g) ■ ■ ■ g) ■ ■ ■ 1 , 

L 

U'l = 1 ■ ■ ■ ■ ■ ■ 1 , (4) 

L 

where Hq and H'j^ are linear operators acting on two dimensional spaces (the 
configuration spaces corresponding to the sites and L, respectively) with bases 
{eo,ei}. 

The number operator in the site i is denoted by n,: 

ni = l(g)---(g)---l(g)n(g)l(g---(g---l. (5) 

i L — i 

n is an operator acting on a two dimensional space with the basis {eo, ei}. The 
matrix form of n in this basis is 

n\ = 5'i5l. (6) 
The evolution equation for the expectation value of an observable Q is 

^^{Q) = {QH), (7) 
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where 

{Q) = SQ^, (8) 

\E' is the (2^+^ dimensional) probabihty vector describing the system and S is 
the covector 

S := s (g) ■ ■ ■ (g) s , (9) 

L 

and 

Sa = 1. (10) 

The system is called autonomous, iff the Hamiltonian is so that the evolution 
of the expectation values of Ui is closed in terms of the expectation values of rij 's. 
In the evolution equation for the expectation value of n^, the expectation values 
of ni_i, Ui, Tii-i-i, {ni^iUi), and (nifii+i) occur. It is seen that the criterion 
that the coefficients of the last two vanish, is 

Sa [H^^f^r® r]°i =0, 

Sa[H,+f,r® r]"' =0, (11) 

respectively, where 

= -5^^ + 51. (12) 

PT|) should hold for all i's, in order that the system be autonomous. So one can 
rewrite it like 

Sa[H^r® r]^'' = 0. (13) 

It is seen that this condition is the same as the corresponding condition for 
uniform lattices, written for each link separately. Provided that this condition 
holds, one arrives at 

^(n-i) = Vi-t^ {ni-i)+9i+fj^ (nj+i) + (Ki_p+!^j+p) (rii) + (Ci-p+o'j+p), < i < L, 

(14) 

where 

ria ■■= Sa {HaY^bt^r^, 

■.= sa iH^y%hr\ 
Sa{H^T\ir\ 

00, 

Sa(i?a)'%0- (15) 
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For the boundary sites (the sites and L), one has 



^(no) = 9^ (m) - + - y^.) (no) + {K\ + cr,,), (16) 



(17) 



3 Relaxation towards the static solution 



The evolution equations governing (?T.i)dy (the difference of (rii) with its static 
value) are the homogeneous parts of (IT4l) . p6)) . and (|17p . which can be written 
like 



^(?^j>dy = ^ ("-i)dy 

The eigenvalues and eigenvectors of the operator h satisfy 



(18) 



Exj =r]j^f_,Xj-i+dj+f_,Xj+i + {Kj^^ + iyj+^)xj, < j < L, (19) 
Exo = 0^,Xl-iH^\+H^\-ly^)xa, (20) 

EXL ^ VL-iiXL-l - {-KL-ii + H'j°i + H']^^o)xL, (21) 



where E is the eigenvalue and x is the eigenvector. The solution to (fTO]) is 



^3 -^j-fJ-y 



where 



and 



Using these, one arrives at 
where 

(PO)) and dlD) lead to 
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Xa — E)ap Xj3, 
E)aP '■= Da-^i ■ ■ ■ £'/3+p- 



Xf, cx 



and 



(22) 
(23) 

(24) 

(25) 
(26) 

(27) 
(28) 
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respectively, where 



Defining 



and 



V 



E 



W := [-rjL-fj^ E - KL-fj^ - VL-f^ + A] , 
it is seen that the eigenvalues can be obtained through 

Things become simpler if one can write Dj as 



D, 



(29) 

(30) 
(31) 
(32) 

(33) 



where is diagonal, and depends on only the parameters corresponding to 
the link a. Putting 

ha b 



(34) 



and 



A, 



Bi 



in (|33l) . exactly similar to [9] one arrives at 

A 



and 



Cla da 



da = 



where 



^a ^a J 
IJa , 
= {-(pka + ■(A) fa, 



(35) 

(36) 
(37) 



(38) 
(39) 
(40) 
(41) 

(42) 
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and (p and ijj are link independent. The parameters k and i> are related to k 
and ly through 

+ '^j+M — '^j-P + '^j+M - ^- (43) 
The consistency of (pS)) to (|^T|) requires that (k — ;>) and [{k + v)"^ — A-qO] be 
constant (link independent). Using 

(k„ + £>„)2 -4r/„6'„ = {na + Vccf - ^Vccda ~ 2 E {k^, + v^) + , (44) 

it turns out that in order that the left hand side be link independent for all of 
the values of E^ both {k + v) and {rj 6) should be link independent. Combining 
these with the link independence of (k — v), it is seen that the necessary and 
sufficient condition for the possibility of the decomposition p3p is that v, 
and (77 0) be link independent. This is one condition more than the conditions 
for superautonomy ( [9]), where it was necessary and sufficient that (k — v) and 
(kv — rjO) be link independent. 

Assuming that ([55)1 holds, one arrives at 



D 



L—fj, /X 



^L~3 fi 3 fi 
i^L — fi 3 fi 



-3 f-L fi ^L — 3 f-i f-i) 



L—fL 3 fi) 



C^L — /.i fi ^L — fi /J,); 



(45) 



where 



dg 



2L 

0/3 ■ 

dp 



Defining Zi through 



^ 'Ha 



it is seen that Zi is link independent and 

da ry 

— — =: Z2 — - 



where 

One then arrives at 



Zi Z2 = sgn(77a 9a). 



(46) 
(47) 

(48) 

(49) 
(50) 



E = ^ + ,,+ y/\^\{Z^' + Z^'), 

where use has been made of the fact that k, and (rjO) are link-independent 
It is seen that (|32|) is equivalent to 

1 ^r-9s 1 



W 



1 



izt' zr') 



1 



\ve\ 



{Z--Z-) 



V^Q, (51) 



6 



or 



= Zi^ — 




(52) 



In the thermodynamic hmit {L — > cx)), ah of the unimodular values of Zi are 
solutions to (|52]). If there are no further solutions, the relaxation time would be 



Tfast 



2Re(V??fi')] 



(53) 



Things change if there are nonunimodular solutions for Z2 as well. In that case, 
one can assume (without loss of generality) that | Z2 1 is larger than one. Then 
([52|) in the thermodynamic limit leads to 



(54) 




or 



Z-. 



c 



A 



(55) 
(56) 



Each of (|55p or ([56|l are acceptable, of course, provided the absolute value of 
the right-hand side is larger than one. Putting these in (f50| . one arrives at the 
following time scales [— Re(i?)]" 



-1. 



-K 
-K 



V + C + C V 



(57) 



One notes that each of t,^ and t\ are larger than rtast, if the right-hand sides of 
([55]) and are larger than one, respectively. So one arrives at the following 
expression for the relaxation time (r) . 



Tfa.t, (-A < ^\) A (-C < 

Tc, (-A < /h^) A (-C > /PI) 

rx, (-A > ^\) A (-C < 

max{Tx,Tc), (-A > /h^) A (-C > /M) 



(58) 



The first case corresponds to the fast phase, where the relaxation time does not 
depend on the boundary rates. The other cases correspond to the slow phase, 
where the relaxation time does depend on the boundary rates. 
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4 Classification of the solvable models 



The local Hamiltonian has 12 independent parameters. Of these 12 parameters, 
only 8 enter the evolution equations of the expectation values of the number 
operators. The criteria for the autonomy of the system are the two equations 
p^ . The system is called superautonomous, if there is a basis corresponding to 
each link, which diagonalizes all one point transfer-matrices corresponding to the 
equation for the time independent configuration of one point functions (hence 
diagonalizing the whole transfer matrix). That condition could be written as 
the existence of matrices A and E such that 

A -S,+mA,E,_^, (59) 

as discussed in [10]. The additional criteria for the system to be superau- 
tonomous and its relaxation be solvable through the transfer-matrix method 
are that k, ly, and {r]9) be link-independent. All of these 5 conditions can be 
written in terms of the following 8 nonnegative independent combinations of the 
rates. 



Aq. 


■■= Sa 


[Ha] 


\a 1 
1 00 


Ba 


■■= Sa 


(Ha. 


\ 1 a 
1 00 


Da 


:= Sa 


(Ha. 


vaO 
) 1 


Fa 


:= Sa 


(Ha. 


\ 1 a 
1 1 


Ga 


Sa 


{Ha, 


vOa 
) 10 


la 


:= Sa 


(Ha. 


\a 1 
) 10 


Ka 


Sa 


(Ha. 


lOa 
) 1 1 


La 


■■= Sa 


[Ha] 


vaO 
) 1 1 



(60) 



The model is autonomous iff 

A-a + Da = + La, 

Ba+Ga^Fa+Ka. (61) 

One has 

Tja — ^a Aa , 

K = —Aa Da, 

V = —Ba — Ga, 
Sa = Fa — Ba- (62) 

Fixing K, Tja, V, and Ga, using (|6T1) and (|62|) . one is still left with two degrees of 
freedom. It should be noted, however, that (ry^ Qa) should be link-independent, 
and (r]a — n), and {9a — v) are nonnegative. Then the parameters entering 
evolution equations are k, u, rja, 9a, Aa, and Ba- The local Hamiltonian, 
however, still contains two more free parameters which do not enter the time 
evolution equations of the expectation values of the number operators. 
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5 Examples 



Consider a lattice each site of which is either empty (o) or full (•). Two examples 
are studied here. 

5.1 The diffused voting model 

The reactions on a link are 

o • — > o o with the rate v — Ca, 

o • — > • o with the rate Cq., 

o • — 5- • • with the rate u — Ca, 

• o — 5- o o with the rate u — ha, 

• o — 5- o • with the rate ha , 

• o — 5- • • with the rate v — ha- (63) 

The rates of injection and extraction of particles in the first (final) sites are a 
and b (a! and 6'), respectively. The rates u and v are link- independent, but the 
diffusion rates, ha and Cq, may be link- dependent. Using and ([^ one 
arrives at 



T] = 




e = 


u, 


K = 


-V, 


V = 


-u. 




0, 


a — 


0, 


c = 


a + b — V, 


A = 


a + b' — u 



(64) 

Although the diffusion rates are link-dependent, all the parameters entering the 
time evolution equation for (ni)dy are link-independent. The time evolution 
equations are exactly the same as those obtained for the corresponding uniform 
model [TT]. There using a different method, it was shown that the dynamical 
phase transition is controlled by the reaction rates at the boundaries. 
Using dnni), dm), and ([Ml), one has 

1 

Tfast 



U V 



T(^— {a + b + u + 



r^^la' + b' + v + —^ ) . (65) 



a + b — V 

u V 
a' + b' ~u 
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It is seen from ([55)) that 




(a + 6) < (v — ^Ju v) 

(a' + h') <{u- Vu^) ■ (66) 
otherwise 

The system is in the fast phase (the last case), when the boundary rates are 
high enough, so that it is the bulk reaction rates that determine the relaxation, 
and goes to the slow phase when the boundary reaction rates are less than some 
critical value. 

5.2 Diffusion and annihilation 

The reactions on a link are 

o • — > o o with the rate w — Bq, 

o • — > • o with the rate Cq, 

• o — > o o with the rate u — ha, 

• o — > o • with the rate ha , 

• • — 5- o • with the rate u — Ca, 

• • — 5- • o with the rate v — ha- (67) 

The rates of injection and extraction of particles in the first (final) sites are a 
and b (a' and b'), respectively. The rates u and v are link-independent, but ha 
and Cq, could be link-dependent, with their product link-independent: 

Cq, ha — m. (68) 

Also each of Ca and ha should be less than or equal to each of u and v, resulting 
in 

m < uv. (69) 

Using p5|) and ((^ one arrives at 



Va = 




9a = 




K — 


-V, 


V = 


—u. 


e = 


0, 




0, 


c = 


a + 


A = 


a' 4 



b — V, 

-b' -u. (70) 
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Using diSl), dST]), and one has 



Tfast = {u + V~2 y/m) \ 



and from ([5 



a + b + u + 



b — V 



Tx - l'^' + b' + v + —^, ) 

\ a' + b' — u J 



(71) 



T = < 



Tfast! {a' + b' > u — -y/m) A (a + 6 > t; 

, (a' + 6' > u - V"^) A (a + 6 < w 

ta , {a' + b' < u - y/m) A [a + b > v 

^max(TA, T^), (a' + 6' > u — -ym) A (a + 6 > ti 



m) 
to) 
to) 
to) 



(72) 



It may be that both (u — \fm) and {v — ^/m) are positive, so unhke the previous 
example it is possible that two phase transitions occur. 



6 Concluding remarks 

Autonomous models are those for which the evolution equations for the ex- 
pectation values of the number operators are closed. Recently the most gen- 
eral autonomous exclusion models with nearest-neighbor interactions and non- 
uniform reaction rates on a one-dimensional lattice were studied. In [10| . using 
a transfer-matrix method, possible static phase transitions of such systems had 
been investigated. Here a similar transfer-matrix method was used to study the 
dynamics and possible dynamical phase transitions of superautonomous mod- 
els. It was seen that superautonomy does not guarantee that the relaxation of 
the system towards its stationary state be solvable using the transfer-matrix 
method. One further constraint is needed, which was found. It was also seen 
that there are two possible phases, regarding the relaxation. In the fast phase, 
the relaxation does not depend on boundary rates. In this phase the relaxation 
time is controlled by only the bulk reaction rates. There may be slow phases as 
well. In these cases, the relaxation does receive contributions from the bound- 
ary rates. A classification was presented for the models solvable through the 
transfer-matrix method, and two examples were explicitly stduied. 

Acknowledgement: This work was partially supported by the research council 
of the Alzahra University. 
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